ApJ accepted [22 February 2013] 

Preprint typeset using I4TgX style emulateapj v. 5/2/1 1 



THE CLUSTER AND FIELD GALAXY AGN FRACTION AT Z = 1 TO L5; EVIDENCE FOR A REVERSAL OF THE 
LOCAL ANTICORRELATION BETWEEN ENVIRONMENT AND AGN FRACTION 

Paul Martini'-^ E. D. Miller,^ M. Brodwin,"* S. A. Stanford,' Anthony H. Gonzalez,* M. Bautz,^ R. C. Hickox,'' 
D. Stern,** P. R. Eisenhardt,* A. Galametz,'' D. Norman,'" B. T. Jannuzi," A. Dey,'" S. Murray,'^ C. Jones,'- & 

M.J.I. Brown'^ 

ApJ accepted [22 February 2013] 

ABSTRACT 

The fraction of cluster galaxies that host luminous AGN is an important probe of AGN fueling processes, 
the cold ISM at the centers of galaxies, and how tightly black holes and galaxies co-evolve. We present a 
new measurement of the AGN fraction in a sample of 13 clusters of galaxies (M > 10'"* Mq) at 1 < z < L5 
selected from the Spitzer/IRAC Shallow Cluster Survey, as well as the field fraction in the immediate vicinity of 
these clusters, and combine these data with measurements from the literature to quantify the relative evolution 
of cluster and field AGN from the present to z ^ 3. We estimate that the cluster AGN fraction at 1 < z < 
1.5 is = 3.0'!^'^% for AGN with a rest-frame, hard X-ray luminosity greater than Lx,h ^ 10"*^ erg/s. This 
fraction is measured relative to all cluster galaxies more luminous than ^(z) + 1 , where ^(z) is the absolute 
magnitude of the break in the galaxy luminosity function at the cluster redshift in the IRAC 3.6/im bandpass. 
The cluster AGN fraction is 30 times greater than the 3(7 upper limit on the value for AGN of similar luminosity 
at z ^ 0.25, as well as more than an order of magnitude greater than the AGN fraction at z ^ 0.75. AGN with 
Lx.H > lO'*^ erg/s exhibit similarly pronounced evolution with redshift. In contrast with the local universe, 
where the luminous AGN fraction is higher in the field than in clusters, the X-ray and MIR-selected AGN 
fractions in the field and clusters are consistent at 1 < z < 15. This is evidence that the cluster AGN population 
has evolved more rapidly than the field population from z ^ 1 .5 to the present. This environment-dependent 
AGN evolution mimics the more rapid evolution of star-forming galaxies in clusters relative to the field. 

Subject headings: galaxies: active - galaxies: clusters: general - galaxies: evolution - X-rays: galaxies - 
X-rays: galaxies: clusters - X-rays: general 



1. INTRODUCTION 

Numerous lines of evidence suggest that there is co- 
evolution, and perhaps a physical connection, between the 
growth of supermassive black holes and the formation of stars 
in galaxies. Perhaps the most striking result is the similar 
rate of evolution of the emissivity from Active Galactic Nuclei 
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(AGN) and star form ation from z ~ 2 to th e present (e.g.lBoyle 
et al.|1998t!F rancesc hinretaL[l 999; Me rloni et al. 2004; SiL 



verman et al.(|2008| ). At the present day, the correlation be- 



tween the masses of supermassive black holes at the centers 
of galaxies and the velocity dispersions o f their spheroid s also 
supports co-evo lution (e.g. |Ferrarese & Merritt 2000, |Geb 



|hardt et aL||2000 ; Tremaine et al. 2002 ) and may inHicate a 
causal connection. Other evidence for a connection between 
black holes and galaxy growth includes that AGN are much 
more comm on in the most luminous starburst galaxies (e.g. 
[Sanders et ~al. 1988; Veilleux et al. 2009) and that even low- 
luminosity AGN are more commonly found in galaxies with 
some young stellar po pulations compared to otherwise similar 
inactive galaxies (e.g. ITerlevich et al. 1 19901 IKauffmann et al.l 

|2DU3t . 

These observational correlations have fueled a lot of inves- 
tigation into the processes that drive matter to accrete onto 
supermassive black holes, as well as form new stars. The 
prevalent, theoretical framework is that the most luminous 
AGN and starburs ts are triggered by major mergers of gas- 
rich galaxies (e.g. [Sanders et al.||1988{ [Barnes & Hemquist] 



1991[ Hopkins et al. 2006 [ Numerous other mechanisms 



have also been proposed to remove angular momentum and 
fuel star formation and black hole growth at lower rates, such 
as large-scale bars, other weakly nonaxisymmetric variations 
in the gravitational potential, m inor mergers, disk instabilities, 
and tu rbulence i n the ISM (e.g.'Simkin et al.'1980''Elmegreen| 
eraLl [l998; Gen zel et al. l 2008; Hopkins & Quataert 201 
ITie observational connection between these mechanisms and 
lower-luminosity AGN is less clear (e.g . |Fuentes-Williams &] 
[Stockel 19881 [Mulchaey & Regan|1997"| [Martini et al.|2003| ), 
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most likely because there are progressively more ways to fuel 
progressively s maller amounts of star formation and black 
hole growth (see Martini|2004| for a review). 

The distribution of AGN in clusters of galaxies relative 
to the field provides some valuable, additional observational 
constraints on fueling processes as a function of luminosity or 
accretion rate, as well as the connection between black hole 
and galaxy growth. This is because additional physical pro- 
cesses impact the availability and transport of the cold gas 
that serves as the primary fuel source for the central black 
hole. These processes include the removal of cold gas via 
ram-pressure stripping (Gunn & Gott 1972), evaporation by 
the hot ISM (Cowie & Songaila 191 1), tida l effects due to the 
cluster potential (Farouki & Shapi ro] 198 lt|M erritt 1983 ) and 
other galaxies (Richstone ,1976t |Moore et al.||1996| ), and gas 
starvation due to the absence of new infall of cold gas ( [Larson | 
|et al."1980). These physical processes have been invoked to 
explain the relative absence of luminous, star-forming galax- 
ies in clusters, the scarcity of substantial reservoirs of cold 
gas, and the large frac t ion of rel atively quiescent, early-type 
galaxies (Gisler 1978 Dressler| [T980; Giovanelli & Haynes] 
[1985) [Dressier et al.|[1999 '). Observational studies of AGN 
in local clusters have similarly found that luminous AGN 
are more rare in cluster galaxies c ompared to fiel d galaxies 
jKauffmann et a l.'2004', l ^opesso & Biviano|2006| ), although 
less-l uminous A GN appea r to be present in compa"able num- 
bers ([Martin i et al. 2002; Miller et al.'2003| [Best et al.|2005a[ 
[Martini et al . 2006; Haggard et al. 2010). 

The different, or at least additional, physical processes that 
influence galaxy evolution in clusters make the cluster envi- 
ronment very well suited to study the co-evolution of super- 
massive black holes and galaxies. This is because the forma- 
tion and evolution of galaxies has proceeded at a different rate 
in clusters relative to the field. For example, the stars in clus- 
ter galaxies appear to have an earlier mean formation epoch 
than field galaxies of simi lar mass (e.g. van Dokkum & Franx 
1996 Kelson et aH1997 ). Star formation in cluster galaxies 



is also observed to mcrease rapidly with redshift (Butcher & 



Oemler'1978';'Saintonge et al.'2008VHaines et al.'2 009(|Hilton 
e^al. 2010 , Tranet al. 2010; Atlee & Martini 2012). Iftheevo- 
lution of AGN in clusters traces the evolution of star-forming 
galaxies in clusters, rather than the evolution of star-forming 
galaxies and AGN in the field, this will be strong evidence 
that AGN and star formation are physically connected, and 
not just a cosmic coincidence. 

The first evidence that substantial numbers of AGN may be 
present in higher-redshift clusters came fro m the discovery of 



three A GN in the z = 0.46 cluster 3C295 by [Dressler & Gunn 
( 1983 ), although subsequent spectroscopic surveys of other 



clust ers at similar redshift s did not find AGN in large num- 
bers ( [Dressier et al.[[T985 ). Later Chandra observations of 
many z > 0.5 clusters did find evidence for AGN through the 
detection of higher surface d ensities of X-ray point sources 
in the fields of these clusters ( Cappelluti et al.|2005[[GiImour[ 



et al. 2009) , although spectroscopi c follow-up observa t ions 
were only obtaine d in a few cases ( Johnson et al.|2003 De- 



marco et al.|2005 ). The first quantitative evidence for a sub- 



stantial increase in the cluste r AGN fraction with redshift was 
presented by ,Eastman et al.[ (2007), who compared the frac- 
tion of spectroscopically confirmed AGN of similar X-ray lu- 
mi nosities in low and high -redshift clusters. A larger study 
by Galametz et al. ( 2009 1 further quantified the increase in 
the AGN traction based on su rface densit y measurements of 
X-ray, MIR, and radio AGN. ,Martini et al.[ ( |2009^ used a spec- 



troscopically confirmed sample to demonstrate that the AGN 
fraction /a increases as (1 for AGN above a hard X- 

ray luminosity of Lx > 10*^ erg s"' hosted by galaxies more 
luminous than M^(z)+ 1, where M|(z) is the absolute magni- 
tude of the knee of the galaxy luminosity function at redshift 
z. This study included a total of 32 clusters from the local 
universe to z ^ 1.3, and included data from many previous 
cluster studies (" Martini et al.|2006[[Eckart et al.|2006 [Martini| 
et al. 2007; Siv akoff et al.|2008] rSeveral more recent studies 
have also id entified AGN in high-redshift clusters and groups 
(Rumbaiigh et al.|2012 [Fassbender et al.|2012[ [Tanaka et al 



ITie rapid rate of AGN evolution is quite similar to the 
evolution of the fraction of st ar-forming ga l axies in clusters 
of fsf oc (1 reported by Haines et al. ( |2009 1, and sug- 

gests the AGN and star-forming galaxy populations evolve at 
similar rates in clusters, although both power-law indices are 
uncertain by approximately ±2. 

The evolut ion of the AGN frac tion in clusters of galaxies 
quantified by [Martini et al.[ ( [2009| ) appears to be substantially 
greater t han the evolution of t he AGN fract ion in the field . 
Work by|Alonso-Herrero et al.l(2008 ) and Bundy et al.|(|2008]) 
demonstrated that the AGN fraction increases by only about 
a factor of two from z ^ 0.5 to z ^ 1.2, which is several 
times smaller than the increase for cluster AGN. This rela- 
tive evolution appears broadly consistent with the behavior 
of sta r-forming galaxies over the same redshift range. Elbaz] 



& [20T2 



et al. (2007j) showed that the fraction of galaxies that are star- 
forming is correlated with local galaxy density at z 1, which 
is a reversal of the anticorrelation observed in the local uni- 
verse. Nevertheless, a direct comparison between field and 
cluster surveys is complicated because they often employ dif- 
ferent selection criteria, such as luminosity in some band, or 
an estimate of the stellar mass, and different AGN luminos- 
ity limits, to establish their host galaxy and AGN samples. 
These selection criteria are important because the AGN frac- 
tion above a given luminosity limit depends on s tellar mass 
(e.g. [Heckman et al.[[2004t [Sivakoff et aL[[2008l lAirdenLl 
I and the evolution of the X-ray luminosity function in- 



dicates that more luminou s AGN were proportionally more 
abundant at higher redshift ( Hasinger et al. 2005 , Barger et aL] 



2005 I, a phenomenon known as "AGN downsizing." 

In this paper we present a new study of a homogenous 
sample of clusters of galaxies at the crucial redshift range of 
1 < z < 1.5 where earlier work implied that the fraction of 
cluster and field AGN would be substantially more similar 
than they are at the present day. These clusters were selected 



from the SpitzerllR AC Shallow Cluster Survey (ISCS, Eisen- 
[hardt et aL[[2008[ ), and have uniform visible, near-infrared, 
and Spitzer observations, deep Chandra observations to iden- 
tify luminous AGN, and substantial photometric and spectro- 
scopic redshift data. We describe these datasets further in the 
next section. We use these data to uniformly select AGN with 
X-ray and MIR criteria, as described in ^ and compare the 
AGN fraction in the clusters with the immediate field envi- 
ronment in ^ to demonstrate the similarity of the field and 
cluster AGN Tractions in this redshift range. In ^we calcu- 
late the cluster AGN fraction for this sample and compare it to 
lower-redshift clusters, then in ^we discuss the relative evo- 
lution of the field and cluster AGN fraction from the present 
day to z - 3. We adopt (VLmMkM = (0.3,0.7,0.7) for the 
cosmological parameters. 

2. DESCRIPTION OF THE DATA 
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TABLE 1 

Properties of the High-Redshift Clusters 



Cluster 


a 


(5 


z. 


ObsID 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


ISCSJ 1429.2+3357 


14:29:15.2 


33:57:08.5 


1.06 


10450 


23 


17.20 


IhCSJ 1432.4+3332 


14:32:29.2 


ii.il.ibM 


1.11 


10452 


34 


17.25 


ISCSJ 1426. 1+3403 


14:26:09.5 


34:03:41.1 


1.14 


10451,7945,6995 


11,41,10 


17.30 


ISCSJ 1426.5+3339 


14:26:30.4 


33:39:33.2 


1.16 


10453 


35 


17.30 


ISCSJ 1434.5+3427 


14:34:30.4 


34:27:12.3 


1.24 


10455 


34 


17.50 


ISCSJ1429.3+3437 


14:29:18.5 


34:37:25.8 


1.26 


10454 


30 


17.50 


ISCSJ1432.6+3436 


14:32:38.4 


34:36:49.0 


1.35 


10456 


32 


17.60 


ISCSJ1425. 3+3428 


14:25:19.3 


34:28:38.2 


1.36 


10458 


36 


17.60 


ISCSJ1433. 8+3325 


14:33:51.1 


33:25:51.1 


1.37 


7946 


40 


17.60 


ISCSJ1434.7+3519 


14:34:46.3 


35:19:33.5 


1.37 


10459 


32 


17.60 


ISCSJ1432.3+3253 


14:32:18.3 


32:53:07.8 


1.40 


10457 


34 


17.65 


ISCSJ1438.1+3414 


14:38:08.7 


34:14:19.2 


1.41 


10461 


lOl(S) 


17.65 


ISCSJ1432.4+3250 


14:32:24.2 


32:50:03.7 


1.49 


10457 


34 


17.75 



Note. — Properties of the clu.sters in the sample. Columns are: (1) Cluster name; (2-3) Right Ascension and Declination (J2000) of the cluster center; (4) Redshift; (5) Chandra 
ObsID of the dataset(s) used in the analysis; (6) Total integration time of the Chandra data; and (7) Vega magnitude at 3.6/im of L* at the cluster redshift from Eisenhardt et al.| 
|2008| . Clusters were observed with the Chandra ACIS-I camera with the exception of ISCSJI438.1+3414 (ACIS-S) 

al., in preparation), where spectroscopic confirmation is de- 
fined to mean that at least five galaxies have redshifts within 
±2000(1 +Zspcc) km s~' of the average spectroscopic redshift 



Our parent cluster sample was sele cted from the IRAC 
Shallow Survey described by Eisenhardt et al. ( 2004| l. This 
survey covers 8.5 de g^ in the NOAO Deep W ide-Field Survey 
in Bootes (NDWFS, [Jannuzi & Deyl[T99 9l) with at least 90s 
of integration time in each of the four IRAC bands. These 
data were supplemented with additional photometry from the 
Spitzer Deep, Wide-Field Survey (Ashby et al. 2009 ), which 
added nine more 30s exposures across the entire area with 
IRAC. Deep Bw, R, and / band data from the Mosaic-1 camera 
on the KPNO 4-m telescope were obtained for the NDWFS. 
Near-infrared images from the FLAMINGOS Extragalactic 
Survey are also available for half of the field (Elston ~et al.| 
|2006). Deeper, near-infrared data have more recently been 
oBtained for the entire field with NEWFIRM (Gonzalez et al. 
in preparation). 

The photometric data were used by Brodwin et al.| ( 2006) l to 
calculate photometric redshifts and redshift probability distri- 
bution functions P{z) with an empirical template-fitting algo- 
rithm. The large AGN and Galaxy Evolution Survey (AGES, 
[Kochanek et al. 2012) of the Bootes field, together with other 
spectroscopic redshift surveys in this region, were used to cre- 
ate training sets and improve the accuracy of the photometric 
redshifts. Based on over 15,000 galaxies with spectroscopic 
redshifts, these photometric redshifts have an uncertainty of 
(7, = 0.06(1 +z) for < z < 1 .5. The photometric redshift cal- 
culations are described in detail in Brodwin et al. (2006). 



and lie within R <2 Mpc of the cluster center. 



Eisenhardt et al. (2008i employed a wavelet analysis tech- 
nique to identify galaxy clusters within the Bootes field for 
the ISCS. The photometric redshift distributions P(z) were 
used to construct weighted galaxy density maps as a func- 
tion of redshift. These density maps in redshift space were 
then convolved with a wavelet kernel and galaxy cluster can- 
didates were identified as peaks in the wavelet-smoothed den- 
sity maps. The significance level for each redshift slice was 
determined with bootstrap resamples of the positions and P(z) 
distributions for the galaxies. Eisenhardt et al. (2008) identi- 
fied a total of 106 cluster candidates at z > 1 with this tech- 
nique and estimated that only ^ 10% may be due to chance 
or projection effects. A number of the se clusters have been 
spectroscopically confirmed to date dStanford et al.||2005 



Our study focuses o n thirteen clusters that we re among the 
most significant from [Eisenhardt et al.| ( |2008| l. Eleven of 
these clusters were targeted as part of a Chandra GO pro- 
gram in 2009 to obtain exposures of 30 ^0 ks (they were also 
included in the lar ger XBootes survey; Murray et al. 2005 
Renter et al.|2005 i. The remaining two clusters have archival 



data with sufficient exposures to inc lude in this study. Two 
other cluster candidates identified by [Eisenhar dt et aL](|2008|l 
also have similar X-ray data, although are not included in this 
analysis. These two clusters are ISCSJ 1427. 9h- 3430, which is 
likely a superposition of ~ 4 groups along the line of sight, 
and ISCSJ 1429. 2H-3425, which appears to be a close pair of 
clusters. We do not include ISCSJ 1427.9H-3430 in our analy- 
sis because it is unlikely to be a massive cluster We do not in- 
clude ISCSJ1429.2H-3425 because the proximity of the cluster 
pair would add uncertainty to our estimate of foreground and 
background contamination. Table [T] lists the basic properties 
of the thirteen clusters, as well as the datasets that we use for 
this study. All of the clusters were observed with the ACIS- 
I camera, with the exception of ISCSJ 1438. Ih-3414, which 
was observed with ACIS-S. The field of view of these ob- 
servations were 16.8' x 16.8' (ACIS-I) or 8.4' x 8.4' (ACIS- 
S). In all cases the cluster was approximately centered in the 
field of view. This field of view is adequate to encompass the 
entire angular extent of these clusters out to approximately 
the r2(K) radius, the radius within which the cluster is a fac- 
tor of 200 overdense relative to the average field value, with 
the exception of some area lost to gaps between the ACIS- 
I chips. These data also include substantial coverage out to 
larger radii, which we use to estimate the field population at 
the cluster redshift. 

We estimated the approximate masses a nd sizes of these 
clusters in several ways. First, Brodwi n" et al.| ( |2007| l 
measured the autocorrelation function of ISCS clusters and 
found ro = 19.14;!^5g h"'Mpc at z = 1, which corresponds 

1014 



to a mean cluster mass of 



M^ 



Based on 



Carl- 



Brodwin et al. 2006; Elston et al.||2()06 ; Eisenhardt et al.| |berg et al.] (1997), we estimate that the typical r2oo raaius 
20081 [Brodwin et al.|2011t |Zeimann et al.|2012| Brodwin et of these clusters is ~ 1 Mpc, which corresponds to 2' 
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at the typical redshift of these clusters. This mean clus- 
ter mass estimate, and thus the inferred size, has bee n con- 
firmed wit h more detailed studies of individual clusters. IBrod-l 
Iwinet al.| (|2011) studied ISCSJ1438.1H-3414 (z= 1.414) and 
ISCSJ1432.4H-3250 (z= 1.487) with deep Chandra observa- 
tions and extensive spectroscopic observations. This study 
showed that ISCSJ1438.1H-3414 has a velocity dispersion of 
747^^Qg km s"' based on 17 members, which corresponds to 
a dynamical mass of log M2oo,rfv/i = 14. S;!;!}^ in solar masses. 
The dynamical mass estimate is in very good agreement with 
the estimate of log M200.X = 14.35;!;q }4 solar masses from 
the cluster's X-ray luminosity. While there is no dynami- 
cal mass estimate for the higher-redshift ISCSJ1432.4H-3250 
cluster, the mass estimated from the X-ray luminosity is sim- 
ilar: log M20[).x = 14.4 ± 0.2 in solar masses. Many of these 



clusters wer e also included in the weak-lensing study by Jee 
and these mass estimates are all above 10'"* M, 



et al. 



(2011 



, 0- 

As we do not have individual mass estimates for all thirteen 
clusters, but do know they are similar based on their selection, 
we adopt a projected radius of 2' for all the clusters. 

There are typically five to fifteen spectroscopically con- 
firmed members in each cluster. Because the spectroscopic 
data are not complete, we supplement these data with clus- 
ter members based on photometric redshifts for some of our 
analysis. We identify these members based on the integral of 
P{z) fromz-0.06(l+z) to z + 0.06(l+z), where z is the clus- 
ter redshift. Galaxies are identified as cluster members if at 
least 30% of the redshift probability distribution is within this 
range and their position is within a projected separation of 2' 
of the cluster center This photometric redshift criterion will 
include foreground and background galaxies. We use a sam- 
ple of galaxies that satisfy the same photometric redshift cri- 
terion, but are projected to lie from 2' up to 10' radius (a phys- 
ical size of r2oo to 5r2oo), to define a field sample and estimate 
the foreground and background contamination. Many studies 
have found evidence that clusters may impact the surround- 
ing field or 'infair region galaxy popu lation at d istances up 
to 5r2oo (e.g. |Patel et al. 2009; Balogh et al.|20091 l. While we 
estimate the field contamination from this region rather than 
the true field, our results below indicate that this distinction is 
unimportant for our analysis. This is likely because galaxies 
that have already been affected by the cluster do not domi- 
nate the surface density of galaxies that satisfy the photomet- 
ric redshift criterion. We also only include galaxies that are 
no fainter than one magnitude below the knee of the luminos- 
ity function, M3 5 + 1 at the cluster redshift. The break in the 
luminosity f unction for cl usters over this redshift range was 
calculated by Mancone et al. (2010). The apparent magnitude 
that corresponds to the break in the luminosity function at the 
redshift of each cluster is listed in Table [T] Figure [T] shows a 
color-magnitude diagram for each cluster, where larger sym- 
bols refer to membership based on spectroscopic redshifts and 
smaller symbols to photometric redshifts. The apparent mag- 
nitudes that correspond to M3 ^ and M3 + 1 for each cluster 
are also shown on each panel. 

3. AGN IDENTIFICATION 

AGN in high-redshift c lusters have pre viously been selected 
based on radio emission pohnson et al.|2 003 ; Galame tz et al. 



2009', 'Gral la et al.|201l|l and X-ray emission (|Johnson et al. 
5003 , East man et al.|2007[[Gaiametz et al.|2009||Martini etaT 



2009 j ). These techniques have unambiguously determined that 
luminous AGN are present in high-redshift clusters because 



only black hole accretion can produce such luminous emis- 
sion at these redshifts. In the first subsection below, we de- 
scribe how we process and then analyze our X-ray observa- 
tions to identify X-ray AGN in clusters and the field and char- 
acterize t heir pro perties. MIR selection based on IRAC color 
selection (|Lacy et al. 2004, Stern et al. 2005, AssefeUiLjSoTO 
Donley et al. 20 lib can also effectively idenmy AGN at these 
redshifts due to tEe very different SED shape of AGN relative 
to normal and star forming galaxies. This is described further 
in the following subsection. 

3.1. X-ray AGN 

The Chandra data were reprocessed with CIAO v4.3 and 
calibration products from CALDB v4.4.1. Point sources were 
identified from a full-resolution image in the observed 0.5-7 
keV energy range using wavdet ect, a wavelet source detec- 
tion tool available in CIAO. We searched at wavelet scales be- 
tween 1 and 32 pixels (0.492 and 15.7 arcsec) to detect a range 
of source sizes and to account for the variable PSF across the 
ACIS field. A wavdet ect threshold of lO"*" was chosen, 
which corresponds to the likelihood of incorrectly detecting a 
source at a given pixel. For the ACIS-I observations, only the 
four ACIS-I detectors were included in the search, therefore 
we expect about four spurious detections in each 2048x2048 
pixel field. One of the clusters (ISCSJ1438.1H-3414) was ob- 
served with ACIS-S centered on the S3 detector; we used the 
same wavdet ect parameters for this observation, only in- 
cluding the single detector Detections of greater than 4 net 
counts were kept in the intermediate source lists, which re- 
sulted in 70 candidate detections of greater than 2-a confi- 
dence in a typical 35 ksec observation. 

One cluster (ISCSJ1426. Ih-3403) has a substantially shorter 
on-axis Chandra exposure (11 ksec). To validate the 
wavdetect results, which were only performed for the 
on-axis observation, we included two additional archival 
datasets, OBSID 7945 (40 ksec) and 6995 (10 ksec). These 
observations were taken with the cluster very close to the 
edge of the field of view, but in opp osite direction s North 
and South. We used acis_extract (Broo s et al.pOlO) , which 
searches for sources in multiple overlapping observations and 
accounts for non-uniform exposures and PSF sizes in the 
searched region. This procedure resulted in no additional 
point source detections within 4' of the cluster center 

We cross correlated our cluster member catalogs and the X- 
ray point source catalogs and identified all cluster members 
with an X-ray source within 2" of the IRAC position. Eleven 
of the cluster galaxies brighter than M3 g + 1 are associated 
with significant X-ray emission. The coordinates, redshifts, 
and apparent magnitudes of these galaxies are listed in Ta- 
ble |2l 

We calculated the flux and luminosity of each X-ray source 
in various observed and rest-frame bands with a F = 1.7 
power-law model, after a correction for Galactic absorption. 
The counts, significance, observed-frame 0.5-8 and 2-8 keV 
fluxes, and rest-frame 0.5-8 and 2-10 keV luminosities are 
provided in Table [3] All of these X-ray sources have rest- 
frame, 2-10 keV luminosities greater than Lxm 

> 10*^ erg/s. 

Only AGN are known to produce point-source emission at 
these luminosities and we consequently classify all these 
galaxies as X-ray AGN. Eight of the 11 X-ray AGN have 
spectroscopic redshifts that confirm they are cluster members, 
including all four with Lx^h > 10'*'* erg/s. The presence of X- 
ray emission was not used as a selection criterion to target 
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Fig. 1 . — Color-magnitude diagrams for cluster members based on spectroscopic (large symbols) or photometric (small symbols) redshifts. All cluster members 
are shown that have photometric uncertainties less than 0.3 mag in the / and 3.6/xm bands and lie within 2' of the cluster center. Error bars are only shown for X- 
ray AGN (blue circles) and MIR AGN (red triangles). The vertical, dotted lines correspond to the apparent magnitudes of Mj ^ andMj g + 1 at the cluster redshift. 
Only cluster members more luminous than MJ^ + 1 are employed in our ana lysis. The vertical, dashed line represents the [3.6] magnitude that coiTesponds to a 
red galaxy at the [4.5] = 17.8 mag limit employed by |Eisenhardt et ar]j2008^ . 
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TABLE 2 

AGN IN THE HIGH-REDSHIFT CLUSTERS 



ID 


Cluster 


Redshift 


a 




S 


/ 


[3.6] 


X/IR 


(1) 


(2) 


(3) 


(4) 




(5) 


(6) 


(7) 


(8) 


J1429I6.1+335537 


ISCSJ1429.2+3357 


1.06" 


14:29:16.1 




+33:55:37.3 


21.30 (0.10) 


15.62 (0.03) 


Both 


J143227.2+333307 


ISCSJ1432.4+3332 


1.15? 


14:32:27.2 




+33:33:07.5 


21.40 (0.08) 


16.01 (0.03) 


X-ray 


J142611.6+340226 


ISCSJ1426. 1+3403 


1.14" 


14:26:11.6 




+34:02:26.0 


21.90 (0.12) 


16.95 (0.03) 


IR 


J142607.6+340309 


ISCSJ1426. 1+3403 


1.12" 


14:26:07.6 




+34:03:09.3 


21.66 (0.09) 


16.44 (0.03) 


IR 


J142602.8+3404()5 


ISCSJ1426. 1+3403 


1.26^ 


14:26:02.8 




+34:04:05.8 


22.30 (0.18) 


16.99 (0.03) 


IR 


J142610.2+340355 


ISCSJ1426. 1+3403 


I.IOP 


14:26:10.2 




+34:03:55.6 


22.33 (0.17) 


17.60 (0.03) 


IR 


J142608.3+340430 


ISCSJ1426. 1+3403 


1.09^ 


14:26:08.3 




+34:04:30.0 


20.83 (0.04) 


16.12(0.03) 


IR 


J142627.5+333912 


ISCSJ1426.5+3339 


1.28^ 


14:26:27.5 




+33:39:12.8 


22.80 (0.27) 


17.08 (0.03) 


IR 


J143430.3+342712 


ISCSJ1434.5+3427 


1.24" 


14:34:30.3 




+34:27:12.0 


21.97 (0.14) 


16.28 (0.03) 


X-ray 


J142915.2+343709 


ISCSJ1429.3+3437 


1.27" 


14:29:15.2 




+34:37:09.2 


21.94 (0.11) 


16.78 (0.03) 


Both 


J143232.7+343538 


ISCSJ1432.6+3436 


1.22^ 


14:32:32.7 




+34:35:38.4 


22.70 (0.28) 


17.36 (0.03) 


IR 


J143237.8+343630 


ISCSJ1432.6+3436 


1.31^ 


14:32:37.8 




+34:36:30.9 


22.46 (0.23) 


16.43 (0.03) 


IR 


J143238.8+343647 


ISCSJ1432.6+3436 


1.28? 


14:32:38.8 




+34:36:47.6 


22.68 (0.26) 


17.20 (0.03) 


IR 


J1425 12.9+342735 


ISCSJ1425.3+3428 


1.21^ 


14:25:12.9 




+34:27:35.2 


21.80 (0.07) 


17.00 (0.03) 


X-ray 


J1425 16.5+342755 


ISCSJ1425.3+3428 


1.36" 


14:25:16.5 




+34:27:55.7 


22.06(0.12) 


17.39 (0.03) 


Both 








1 zt-9S-9n ^ 

IM-.Z J .ZU. J 










Both 


J142520.2+343014 


ISCSJ1425. 3+3428 


1.36^ 


14:25:20.2 




+34:30:14.3 


23.48 (0.29) 


17.31 (0.03) 


IR 


J1425 19.8+343024 


ISCSJ1425.3+3428 


1.17^ 


14:25:19.8 




+34:30:24.8 


22.23 (0.09) 


16.90 (0.03) 


IR 


J143351.5+332645 


ISCSJ1433. 8+3325 


1.37" 


14:33:51.5 




+33:26:45.8 


20.19(0.03) 


16.12(0.03) 


Both 


J143445.7+351921 


ISCSJ1434.7+3519 


1.37" 


14:34:45.7 




+35:19:21.7 


22.69 (0.18) 


16.94 (0.03) 


IR 


J143450.0+351958 


ISCSJ1434.7+3519 


1.43P 


14:34:50.0 




+35:19:58.8 


22.94 (0.28) 


17.34 (0.03) 


IR 


J143216.4+325434 


ISCSJ1432.3+3253 


1.39" 


14:32:16.4 




+32:54:34.1 


20.43 (0.05) 


16.35 (0.03) 


Both 


J143217.1+325235 


ISCSJ1432.3+3253 


1.27^ 


14:32:17.1 




+32:52:35.1 


22.88 (0.28) 


17.49 (0.03) 


IR 


J143218.1+325315 


ISCSJ1432.3+3253 


1.40" 


14:32:18.1 




+32:53:15.9 


21.85 (0.16) 


16.69 (0.03) 


IR 


J143816.8+341440 


ISCSJ1438.1+3414 


1.41" 


14:38:16.8 




-1-34:14:40.3 


22.41 (0.16) 


17.05 (0.03) 


Both 


J143817.4+341337 


ISCSJ1438.1+3414 


1.46'' 


14:38:17.4 




+34:13:37.6 


22.52 (0.16) 


17.58 (0.03) 


IR 


J143802.7+341548 


ISCSJ1438.1+3414 


1.42^ 


14:38:02.7 




+34:15:48.2 


21.73 (0.23) 


16.61 (0.03) 


IR 


J143217.1+325055 


ISCSJ1432.4+3250 


1.58^ 


14:32:17.1 




+32:50:55.1 


22.60 (0.19) 


17.49 (0.03) 


Both 


J143224.8+325005 


ISCSJ1432.4+3250 


1.31'' 


14:32:24.8 




+32:50:05.0 


22.59 (0.14) 


17.08 (0.03) 


IR 


Note. — AGN in 


he high-redshift clusters. Columns 


are: (1) AGN ID; (2) Cluster name; (3) AGN Redshift, either photometric (p) or spectroscopic (s); (4—5) Right Ascension 


and Declination (J2000); (6-7) / and [3.6] (Vega) mag; and (8) AGN selection via X-ray, IR, or Both criteria. For AGN with photometric redshifts, the redshift listed in column 3 


may not agree with the cluster redshift because it corresponds to the peak of the redshift probability distribution function, while membership is based on the fraction of the integrated 


probability at the cluster redshift. See ^for further details. 
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X-RAY Properties of the Cluster AGN 








ID 


Cluster 


Redshift 


Counts Si 


g 


Fo.s-s 




^0.5-8 


-^-2-10 


(1) 


(2) 


(3) 


(4) (5) 


(6) 


(7) 


(8) 


(9) 


J142916.1+335537 


ISCSJ1429.2+3357 


1.06" 


89 44 





6.71 (0.70) 


4.04 (0.42) 


3.20 (0.33) 


2.44 (0.26) 


J143227.2+333307 


ISCSJ1432.4+3332 


1.15'' 


6 3 


5 


0.25 (0.09) 


0.15 (0.06) 


0.14(0.05) 


0.11 (0.04) 


J143430.3+342712 


ISCSJ1434.5+3427 


1.24" 


11 6 





0.49 (0.13) 


0.29 (0.08) 


0.33 (0.09) 


0.26 (0.07) 


J142915.2+343709 


ISCSJ1429.3+3437 


1.27" 


75 37 





4.11 (0.47) 


2.48 (0.28) 


2.91 (0.33) 


2.28 (0.26) 


J1425 12.9+342735 


ISCSJ1425.3+3428 


1.21'' 


4 2 


4 


0.20 (0.08) 


0.12(0.05) 


0.12(0.05) 


0.10(0.04) 


J1425 16.5+342755 


ISCS J 1425. 3+3428 


1.36" 


25 12 


8 


0.99 (0.18) 


0.59(0.11) 


0.81 (0.15) 


0.65 (0.12) 


J142520.3+342942 


ISCSJ1425.3+3428 


1.36" 


25 12 


4 


0.90 (0.17) 


0.54 (0.10) 


0.72 (0.17) 


0.58 (0.14) 


J143351.5+332645 


ISCSJ1433. 8+3325 


1.37" 


81 39 


1 


2.65 (0.28) 


1.60(0.17) 


2.21 (0.24) 


1.76 (0.19) 


J143216.4+325434 


ISCSJ1432.3+3253 


1.39" 


63 25 


9 


2.55 (0.31) 


1.54 (0.19) 


2.21 (0.27) 


1.75 (0.21) 


J143816.8+341440 


ISCSJ1438.1+3414 


1.41" 


10 4 


5 


0.30 (0.05) 


0.18(0.03) 


0.26 (0.04) 


0.21 (0.03) 


J143217.1+325055 


ISCSJ1432.4+3250 


1.58'' 


9 5 





0.43 (0.14) 


0.26 (0.08) 


0.49(0.15) 


0.40 (0.13) 



Note. — X-ray properties of the AGN in the high-redshift clusters. Columns aie: (1) AGN ID; (2) Cluster name; (3) AGN Redshift, either photometric (p) or spectroscopic (s); 
(4) Counts; (5) Significance of the X-ray detection; (6-7) Unabsorbed flux in the observed- frame 0.5-8 keV and 2-B keV bands in units of 10"^"^ erg cm~^ s~'; (8-9) Unabsorbed 
luminosity in the rest-frame 0.5-8 keV and 2-10 keV bands in units of 10'^'^ erg s~^ . 
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candidate cluster members for spectroscopy; however, two of 
the spectroscopically confirmed AGN were selected based on 
their MIR colors, as described in the next subsection. 

Our X-ray data do not have uniform depth over the area 
subtended by each cluster and field region. The nonunifor- 
mity is due to the gaps between the ACIS-I detectors and how 
well the clusters were centered on the detectors. We used ex- 
posure maps for each cluster to quantify the fraction of the 
cluster and field area with sufficient depth to identify AGN. 
We identified a sensitivity threshold for each exposure map 
that corresponds to an AGN with a rest-frame 2-10 keV lumi- 
nosity of 10"*'* erg/s at the cluster redshift. The fraction of the 
cluster and field regions above this threshold for each cluster 
is used to quantify the X-ray AGN surface density and X-ray 
AGN fraction in the following sections. 

3.2. MIR AGN 

We used our IRAC data to identify candidate cluste r AGN 
with the MIR color-selection criteria defined by Stern et al.| 
(12005 ) and refer to these as MIR AGN. Color-color diagrams 
tor the clusters are shown in Figure |2] Galaxies were only 
classified if the photometric uncertainties were less than 0.3 
mag in all four IRAC bands. There are 27 MIR AGN, of 
which 12 have spectroscopic redshifts. Eight of the 11 X- 
ray AGN were also identified as MIR AGN, including all 
the X-ray AGN more luminous than Lx.h > 10^*4 erg/s. The 
incomplete overlap between the MIR AGN and the X-ray 
AGN is consistent with previous studies, which have found 
similar results for the field (Hick ox et al.||2009] ) and in low- 
redshift clusters ( |Atlee et al.||201 l| l. Unlike the case for 
X-ray AGN, J142916.1H-335537 in ISCSJ1429.2H-3357 and 
J143816.8H-341440 in ISCSJ1438.1H-3414 were selected for 
spectroscopic observations as cluster MIR AGN candidates. 
Both of these MIR AGN are also X-ray AGN. 

4. RADIAL DISTRIBUTIONS 

The projected distribution of all galaxies consistent with 
each cluster's redshift, including the X-ray and MIR AGN, 
is shown in Figure |3] The adopted angular size for each of the 
clusters is 2', which is approximately the value of r2()o (~ 1 
Mpc) at the redshift of these clusters (see Field galax- 
ies consistent with the cluster redshift are shown within an 
8' X 8' box centered on each cluster to illustrate the environ- 
ment in the immediate vicinity of these clusters. We use these 
galaxies, including AGN, from /? = 2' — > 10' to measure the 
surface density of field galaxies and quantify the amount of 
foreground and background contamination in the cluster sam- 
ple. A fraction of the cluster and field area does not have com- 
plete X-ray coverage due to chip gaps and the exact placement 
of the cluster center relative to the center of the ACIS field of 
view. We corre ct for the non-uniform X-ray coverage as de- 
scribed in §3.1| 

As luminous AGN are too rare to study their distribution 
within individual clusters, we stack the cluster catalogs to 
measure their radial distribution. The total number of sources 
per 0.5' bin as a function of distance from the cluster center is 
shown in Figure]?] The solid lines correspond to all galaxies 
that are at the cluster redshift. The dotted lines only include 
galaxies with spectroscopic redshifts. The almost complete 
lack of galaxies with spectroscopy outside of 2' is because 
these galaxies were generally not spectroscopic targets. Be- 
cause the spectroscopic coverage is fairly high in the clusters, 
we only use AGN with spectroscopic redshifts to compare to 



AGN in clusters at other redshifts. We use photometric red- 
shifts to estimate the total number of cluster galaxies, as well 
as to estimate the relative AGN fractions of X-ray and MIR 
AGN in clusters and the field. 

Figure]5]shows the surface density of all galaxies {left), and 
X-ray andMIR AGN (right). The surface densities for all 
galaxies and for MIR AGN are simply the raw counts shown 
in Figure ]4] divided by the total area of each annulus. The 
surface density of X-ray AGN is computed from the total area 
in each radial bin that is above the fixed luminosity sensitivity 
threshold described in 33] The surface density of galaxies, X- 
ray AGN, and MIR AGN all asymptote to constant values by 
approximately 2' from the cluster center, which is consistent 
with the adopted radius of 2' for these clusters. We estimate 
that small uncertainties in our choice of 2' for the radius of 
these clusters will not impact our results. 

We use the data from 2' — > 10' to estimate the field sur- 
face density of all galaxies and MIR AGN, while we use the 
data from 2' — > 6' for all galaxies that have X-ray coverage 
and X-ray AGN. The radial range for the X-ray coverage is 
smaller because of the size of the X-ray images and the de- 
terioration of the point spread function further off axis. The 
surface densities for all galaxies, all galaxies with X-ray cov- 
erage. X-ray AGN, and MIR AGN are 0.90, 0.99, 0.04, and 
0.07 arcmin"-, respectively. We use the field density to cal- 
culate the foreground and background contamination within 
the clusters, and then calculate and subtract an estimate of the 
contamination from the total number of galaxies within a pro- 
jected radius of R = 2' to estimate the total number of cluster 
members. X-ray AGN, and MIR AGN in the cluster sample. 
In the case of X-ray AGN, this includes a correction that ac- 
counts for the fact that not all the projected cluster and field 
area is above the X-ray luminosity threshold. We estimate that 
there are a total 150 galaxies in these thirteen clusters above 
the galaxy luminosity threshold and 136 of these galaxies lie 
within the area above our X-ray luminosity threshold. This 
surface density extrapolation also leads to an expectation of 
five X-ray AGN and 15 MIR AGN in these clusters. These 
values are consistent with the eight and 12 spectroscopically 
confirmed X-ray and MIR AGN in these clusters (eight AGN 
are common to both samples) listed in Tables]2]and]3] 

The fractional surface densities of X-ray and MIR AGN in 
the cluster and field samples are consistent with one another. 
The X-ray and MIR AGN fractions in the clusters and the 
field are shown in Figure ]6] Note that these fractions do not 
correspond to all AGN above a fixed luminosity threshold, 
but are instead calculated from all AGN above our detection 
threshold in galaxies more luminous than g + 1 ■ An alter- 
native way to express this is to note that the cluster galaxies, 
cluster galaxies with X-ray coverage. X-ray AGN, and MIR 
AGN correspond to similar overdensities of (S- (S))/ (S) = 
1.02,0.87,0.82, and 1.25, respectively. We discuss the signif- 
icance of the similar cluster and field AGN fractions in ^ 

There are three biases that influence our estimates. First, 
the cluster sample has been identified via spectroscopic red- 
shifts, or photometric redshifts if no spectroscopic data are 
available, while the field sample has almost exclusively been 
identified via photometric redshifts. As a result of the spec- 
troscopic observations, the surface density of foreground and 
background contamination within the clusters is lower than it 
is outside of them. We have compared the number of cluster 
members within 2' based on photometric redshifts alone to 
the number when spectroscopic data are included and find the 
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Fig. 2. — IRAC color-color plots for members of the clusters listed in Table^(o// points). Error bars are only shown for X-ray AGN (blue circles) and MIR 
AGN (red triangles). Large symbols correspond to spectroscopically confirmedcluster members. Small symbols correspond to photometric redshifts. Galaxies 
are only shown if they are projected to lie within 2' of the cluster center, have photometric uncertainties less than 0.3 mag in all four IRAC bands, and they are 
brighter than M*g(z) + 1, where M^^(z) is the [3.6] magnitude of the break in the luminosity function at the cluster redshift. 
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Fig. 3. — Positions of all cluster galaxies out to a projected distance of 2' r200. dashed circle) and field galaxies consistent with the cluster redshift within 
an 8' X 8' box centered on the cluster. Only galaxies with photometric uncertainties less than 0.3 mag in the / and 3.6/^m bands are shown. Symbols are as in 
Figures[T]and[2] 
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Fig. 4. — Stacked histograms of the galaxies shown in Figureplas a function of clustercentric distance out to 4'. (Left) Stack of all galaxies with photometric 
or spectroscopic redshifts (solid line) and only spectroscopic reoshifts (dotted line) consistent with cluster membership as a function of projected clustercentric 
distance. (Right) As at left for the X-ray AGN (blue lines) and MIR AGN (red lines) as a function of projected clustercentric distance. 
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Fig. 5. — Same as Figure[4] except out to 10' from the cluster centers and normalized by the ai'ea of each annulus to show the density of sources per square 
arcminute. The dashed, horizontal lines correspond to the median field surface densities calculated from 2' to 10' (for all galaxies and IR AGN), or from 2' to 6' 
(for X-ray AGN) from the cluster center. 

change is ^ 10%, that is we may have slightly underestimated (see ^3.2[ ) and both are also X-ray AGN. These constitute only 
the total number of cluster galaxies that form the denominator two out of 1 1 X-ray AGN and two out of 27 MIR AGN, so 
of the cluster AGN fraction (this will not affect the relative this selection has a minor impact on the relative cluster and 
AGN fraction in the field and clusters). The second bias is field AGN fractions. This bias does not impact the numera- 
that two MIR AGN were selected as targets for spectroscopy tor of the AGN fractions that are derived from spectroscopic 
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Fig . 6 . — Relative fraction of X-ray AGN (blue) and MIR AGN (red) in the 
field and cluster samples. The en'or bars correspond to the 90% confidence 
limits calculated from the relative Poisson uncertainty in the number of AGN 
in each subsample. The field sample is comprised of AGN located 2' — 10' 
(IR AGN) and 2' - 6' (X-ray AGN) from each cluster that have been selected 
in an otherwise identical manner to the cluster sample. 

redshifts alone. The third bias is that our field region only ex- 
tends to a projected separation of 3r2()o (X-ray AGN) to 5r2oo 
(IR AGN). The region out to 5r2oo is sometimes referred to 
as the "infall region" as galaxies at these distances may have 
been "pre-processed" by membership in an inf alling group or 
have even a lready passed through the cluster (Diaferio et aL 
IDoniPateTet al. 2009; Balogh et al. 2009; Bahe et al. 201 



and thus may not be representative of the true field population. 
We expect most of the field sample is representative of the true 
field because the photometric redshifts will include many true 
foreground and background galaxies in the field sample; nev- 
ertheless, our field estimate may be somewhat biased by the 
infall population. 

5. EVOLUTION OF THE AGN FRACTION IN CLUSTERS OF 
GALAXIES 

Eight of the 1 1 X-ray AGN in these clusters have spectro- 
scopic redshifts that confirm their cluster membership. All 
of these X-ray AGN have rest-frame, hard X-ray luminosi- 
ties greater than Lx.h > 10^^ erg/s and four (all with spectro- 
scopic redshifts) have Lx.h > 10'*'* erg/s. Our X-ray data are 
sufficiently sensitive that we should have detected all AGN 
with Lx.h > 10'*'* erg/s at the redshifts of these clusters, with 
the possible exception of AGN in the chip gaps of the ACIS 
camera. We have examined exposure maps of these fields and 
conclude that chip gaps affect at most five of the 350 galaxies 
that are within the area of these clusters and have photometric 
or spectroscopic redshifts consistent with cluster membership. 
Therefore a negligible fraction of X-ray AGN are missed in 
the chip gaps. In |4j we estimated that there are a total of 136 
cluster members in galaxies more luminous than M3 + 1 with 
complete X-ray coverage (after we correct for foreground and 
background contamination). We therefore estimate that the 
AGN fraction is /a(Lx,h > 10"*"*) = 0.030![|:[]?t The uncertain- 
ties are Poisson errors that cor respond to la Gaussian confi- 
dence intervals |Gehrels|1986[ ). 

We are only sensitive to sources as dim as Lx.h ^ 10'*^ erg/s 
in a small minority of the clusters. As the luminosity function 
of cluster AGN is not well known at these luminosities, and 
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Fig. 7. — Evolution of the X-ray AGN fraction in clusters from z = to 
z = 1.5 for hard X-ray luminosity thresholds of Lx.h > erg/s (red) and 
Lx.h > 10*'' erg/s (blue). The error bars correspond to Icr Gaussian errors 
and the upper limit at z ^ 0.25 corresponds to a 3a upper limit. The lower 
limit at z ~ 1.25 is due to incompleteness. 

variations in the spectral energy distributions and intrinsic ab- 
sorption are similarly not well known, we do not attempt to 
model and correct these data to estimate the total number of 
X-ray AGN with Lx.h > 10'*l histead we use our data to es- 
tablish a lower limit on the cluster AGN fraction for sources 
withLx,H > lO'*^ erg/s. There are eight spectroscopically con- 
firmed AGN in this category (and eleven total), for a lower 
hmit offAiLxM > 10"*^) > 0.059!j;:|J^?(0.082!;]:|;^2). The AGN 
fractions based on the spectroscopically confirmed AGN are 
shown in Figure |7] near z = 1 .26, the median redshift of the 
cluster sample (although offset slightly in redshift for clarity). 

Both of these measurements imply that luminous X-ray 
AGN are more common in clusters at z > 1 compared to 
lower redshifts, as had been indicated by previous work at 
Z < 1 (Eastman et al.||2007| Galametz et al. 2009, Ma£] 
|tini et al.,.2009j . [Martini et al.. ( ,2009) calculated the AGN 
fraction and evolution in a similar manner to this study for 
32 clusters at 0.05 < z < 1.27 (only three at z > 1) and 
Lx.H ^ 10"*^ erg/s. They found that the AGN fraction was 
fA(Lx,H > lO'*^) = 0.00134 for 0.05 < z < 0.4 and 0.0100 
for 0.4 < z < 1.27, where the uncertainty in the measure- 
ment at z < 0.4 is dominated by the presence of only two 
AGN in the 17 clusters in this redshift range. We have com- 
bined the 19 clusters at z < 0.5 presented in [Martini et al.| 
(2009) with the more recent study by Haines et al. (20O]l 
of AGN in clusters at 0.16 < z < 0.29 to construct a com- 
bined sample of 44 cl usters and 3869 + 27 02 = 6571 cluster 



Martini et al. 



galaxies at z < 

(two with 0.4 < z < 0.5 ) with Lx.h > 
with Lx.h > lO'*'* erg/s). 



( [2009] found four AGN 
10^^ erg/s (and none 



Haines et al. 



(20121 quote X-ray lu- 



minosities in the 0.3-7 keV band, rather than 2-10 keV. We 
used the ratio of the 0.3-7 keV band to the 2-10 keV band 
for a r = 1.7 power law to scale their results to 2-10 keV. 
This yields three AGN in their sample with Lx,h > lO'*^ erg/s 
and none above lO'*'* erg/s. The combined sample has AGN 
fractions of MLx,h > 10«) = 0.00107!»;J]0»^7 and /^(Lx,^ > 



Many of the clusters at 0.2 < z < 0.4 studied by Hart et al. i2009| over- 
lap with either the Martini et al. ( 2009 1 or Haines et al. ( 20 1 2 1 sam ple. Only 
one clu ster (A1689) overlaps between |Martini et al.n2009) and |Haines et al.| 
(20T2l. 
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10"*^) < 0.00101 (a 3(7 upper limit). This point and upper 
limit are shown in Figur e [7| near z = 0.25. At intermediate 
redshifts of 0.5 < z < l, |Martini et al.| ( |2009] l found 13 AGN 
with Lx,H > lO'*'' erg/s, including two with Lx.h > 10"^^ erg/s, 
in ten clusters and an estimated population of 1734 galax- 
ies. The AGN fractions are fAiLx.H > ^0^^) = 0.0075!;j:|j;]^^ 
and fA(LxH > 10*"^) = 0.0012!°;°°'^. These points are shown 
in Figure j7] near z = 0.75. Note the luminosity thresholds 
employed to define the galaxy samples in these studies are 
similar to ours: Martini et al. (2009) adopted a threshold of 
M^(z) + 1 and |Haines et al (2012) adopted M|. + 1 .5. 

The evolution in the AGN fraction corresponds to a fac- 
tor of at least 45 for Lx.h > 10"*^ erg/s AGN from z ^ 0.25 
to z ^ 1.25 and a factor of at least 30 for Lx,h > 10"*"^ erg/s 
AGN. In the former case, the increase is a lower limit due 
to incompleteness at z > 1 . In the latter case the increase is 
a lower limit due to the absence of any AGN in this lumi- 
nosity range in the 44 clusters that comprise the low-redshift 
sample. For the two luminosity and redshift ranges without 
limits on the fractions, the evolution in the AGN fraction is 
also pronounced. For Lx,h > lO'*^ erg/s AGN, the fraction in- 
creases by a factor of seven over the range 0.25 < z < 0.75. 
For LxH> 10"^'* erg/s AGN, the fraction increases by a factor 
of 25 over the range 0.75 < z <C 1 .25. This substantial increase 
is in good agree ment with the i ncrease ove r this same redshif t 
range found by Galametz et al.| ( ,2009) and Hart et al. ( 201 l| l, 
although for a somewhat different range in luminosities. 

Two complications in the interpretation of the evolution of 
the AGN fraction are the extent to which the high-redshift 
clusters are equivalent to the progenitors of the low-redshift 
clusters and the incompleteness and evolution of the cluster 
galaxies. If the AGN fraction is a strong function of clus- 
ter mass, in addition to redshift, then mass dependence could 
manifest as redshift depende nce. Based on previous analysis 
of the 1 < z < 1 .5 sample ( Brodwin et al 
these clusters are expected to be 
progenitors of present-day 10'^ Mq clusters such as Coma. 



S revious analysis 
2007l[20TT| |Jeel 
e the high-redshift 



The cluster sam ples in the range < z < 1 studied by Mar- 
tini et al. ( 2009 1 are similarly massive and consistent with the 



same population of clusters in th e local univ erse. Similarity 
to Coma was explicitly used by [Hart et al .l(2009 ) to select 
their sample, which overlaps with the other low-redshift clus- 
ter samples. To the extent that there is a bias in these data, the 
higher-redshift clusters tend to be slightly more massive. As 
there is an anticorrelation be tween luminous AGN a nd envi- 
ronment in the local universe ( |Kauffmann et al.|2004[ ), this an- 
ticorrelation would produce an underestimate of the true rate 
of evolution. 

The other complication originates in the comparison of the 
cluster galaxy populations across different surveys and from 
the local universe to z = 1.5. In this study we have focused on 
the AGN fraction for cluster galaxies more luminous than one 
magnitude below the knee of the luminosity function at the 
cluster red shift. This c hoice was motivated by the previous 
study by Martini et al.| ([2009), who adopted the assumption 
that Ml(z) = Mg(0)-z. [TTaines et al.|(|2()12|l adopted a simi- 



lar threshold for the galaxy population, although their galaxy 
luminosity threshold was 1.5 magnitudes below the knee of 
the luminosity function (defined in the /T-band: M|.+ 1.5). 
We estimate that the half magnitude difference for the galaxy 
population is a minor effect compared to the Poisson uncer- 
tainties, as the z < 0.5 AGN fractions measured by these two 
studies are consistent. At intermediate redshifts 0.5 < z < 1 



the cluster spectroscopic data are incomplete for galaxies (al- 
though not for the much smaller number of X-ray sources), 
and there is thus a substantial correc tion to estimate the t otal 
galaxy population of some clusters. |Martini et al.| ( |20()9| in- 
vestigated this effect and estimated that it introduced on order 
a factor of two uncertainty into the AGN fraction estimate at 
z > 0.5, which is much smaller than the observed evolution. 
While the present sample at 1 < z < 1 .5 similarly suffers from 
spectroscopic incompleteness, the high-quality photometric 
redshifts, substantial sample of field galaxies, and availabil- 
ity of some spectroscopic follow-up substantially mitigate this 
uncertainty such that the cluster galaxy population is proba- 
bly more reliably known for the present sample than for the 
sample at 0.5 < z < 1. 

6. DISCUSSION 

The field X-ray AGN fraction also increases over the same 
redshift range of < z < 1.5 where the cluster AGN fraction 
has in creased by factors of at least 30-45 (e.g. Ueda et al. 
2003| l. At z 1.25, we showed in ^and Figure |6] that the 
cluster and field AGN fractions were comparable and both 
~ 3-4% for AGN with hard X-ray luminosities greater than 
a few times 10'*"' erg/s. Because our field measurement is from 
the immediate vicinity of massive clusters, some of the field 
galaxies may have already been process ed through the clus- 
ter (e.g. IDiaferio et a l][2dni] [Patel et al. 2009, BaloghenT 
|2009 i Bahe et al. 20l5jl and the field sample may not be repre- 
sentative of the true field population. Another estimat e of the 
field A GN fraction at these redshifts was presented by |Bundy] 
'et al.' (2008) from an analysis of galaxies and AGN in the 
DEEP2 survey. They found that the AGN fraction is - 1 - 3% 
for AGN with Lx^h > 10"*^ erg/s in host galaxies with stellar 
masses of ~ lO^*'^ at 1 < z < 1.4. This fraction appears 
consistent with our measurement for the field fraction around 
the clusters, although our data do not extend to Lx.h = lO'*'' 
erg/s for most clust ers and the stellar mass range quoted by 
[Bundy et al. ( |2008] l is not an exact match to our luminosity 
threshold at 3.6/im. 

In the local universe, very luminous AGN are rarely found 
in the field and very r arely found in clusters. An SDSS study 
by Kauffma nrTeTaD ( 2004) found that AGN with L[0 III]> 
10' L0 were approximately three times rarer in dense envi- 
ronments relative to less dense environments, where this lu- 
minosity threshold is approximately comparable to an AGN 
with Lx.h = 10^'' erg/s. A factor of three to four decrease 
in the AGN f raction with the same luminosity threshold was 
also found by Best et al. ( 2005b ) over an order of magnitude 
change in local density with data from SDSS. However, both 
studies also find little environmental dependence in the AGN 
fraction when only lower-luminosity AGN were examined. 
Note that for both of these studies, the high density regions 
are still not as dense as the centers of rich clusters of galax- 
ies, which are embedded in very rare and very massive dark 
matter halos (on order lO''' Mq) in the local universe. 

The local field X-ray AGN fraction was measured by|Hag-| 
Igard et aT] ( |2010[ ) based on data from the Cha ndra Multi- 
wavelength Project (ChaMP; Green et al.| [200?) and SDSS, 
although they do not use hard X-ray luminosities and probe 
a somewhat different range in galaxy luminosity and redshift. 



Haggard et al. ( 2010 1 compared their re sults to previous clus- 
ter studies ( [Martini et aL][2006l [2007) 1 and found the field 
and cluster AGN fractions are the same for low-luminosity 
AGN {Lx > 10"*'"^^ erg/s). Other studies have reached simi- 
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Fig. 8. — Evolution of the X-ray AGN fraction in clusters {solid symbols) and the field (open symbols) from z = to z ~ 3. All of the cluster measurements 
have been reproduced from Figure^ The field AGN fractions at z ~ 0.3 for Lx,h > 10*' erg/s (open, red circles) and Lx,h > 10'*'' erg/s (open, blue circles) are 
from D. Haggard (private communication, see also Haggard et al.pOlO) . The formal uncertainties on these two field fractions are smaller than the size of the 
circles. The field AGN point at z ~ 1.25 (open blue circle) is scale d from the clust er point by the ratio of the field and cluster fractions shown in Figurel6] The 
z = 2.3 protocluster and neighboring field fractions from Digby-No rth et al^j2010^ (filled and open green hexagons, respectively) and the z = 3.09 protocluster 
and neighboring field fractions from |Lehmer et ah]|2009^ (filled and open green pentagons, respectively) are also shown. These points are described further in 



lar conclusions, namely the fraction of low-luminosity AGN 
is comparabl e in clusters , groups, and the field in the lo- 
cal universe ( |Sivakoff et a l. 2008 ; Arnold et al.|2009| |Miller| 
et al. 2012| l. The sample studied by Haggard et al. (2010 1 
also includes some higher-luminosity AGN that are more di- 
rectly comparable to the AGN considered here. D. Hag- 
gard {private communication) has computed the field AGN 
fraction for similar X-ray and galaxy luminosity limits for 
< z < 0.6 and found /a(Lx,h > 10"*^) = 0.0064!;;{;|;[!^ and 
fAiJ-x.H > 10'*'*) = 0.001l!j]:|]jS|]^. These points are shown in 
Figure [8] (note the formal uncertainties are smaller than the 
points) and demonstrate that the Lx.h > 10"*^ erg s"' field 
AGN fraction is six times higher than the cluster value and 
the Lx.h > lO'*'* erg s"' field AGN fraction is consistent with 
the 3(7 upper limit for the cluster AGN fraction. While the 
two field fractions at z ^ 0.3 are calculated with the ChaMP 
survey's definition of the hard band of 2-8 keV, rather than the 
2-10 keV band adopted throughout the rest of this paper, this 
is a very minor difference. The relative field and cluster AGN 
fractions in the local universe and at z ^ 1 .25 show that while 
luminous AGN are anticorrelated with local density in the lo- 



cal universe, this is no longer the case at z ^ 1 .25. The masses 
of these high-redshift clusters are also in the range expected 
for the progenitors of the local clusters. 

Studies at even higher redshift support this trend and sug- 
gest that the present-day anticorrelation has reversed by z > 2. 
Chandra observations of three protoclusters at z > 2 have 
revealed luminous AGN associated wit h PKS 1138-262 at 
z = 2.16 ( |Pentericci et al.|2002HCroftera l. 2005 ), the z = 2.3 



protocluster in the field of QSO HS 1700H-643 (Digby-North 
et al. 2010), and th e z = 3.0 9 SSA22 protocluster (Lehmer 
et al. 2009). Pentericci et al. (2002 ) compared the number of 
sources toward PKS 1 138-262 and found an excess of ~ 50% 
compare d to expectations from the AGN space density at this 
redshift. [Lehmer et aT] ([2009^ and i Pigby-North et al.| ( |2010) l 
both measure the AGN fractions o f the protoclus t ers an d in 
field samples at the same redshift. |Lehmer et al. (2009) de- 
tected X-ray emission from six LBGs and five LAEs (ten 
unique sources) toward the SSA22 protocluster at z = 3.09 
with a 400ks Chandra observation. These sources have X-ray 
luminosities of 3-50 x lO'*^ erg/s in the rest-frame 8-32 keV 
band. They measure AGN fractions of 9.5!^-i^% and 5.l!f j% 
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in LBGs and LAEs, respectively. These AGN fractions are ( [Dressier et al.|1999l l. 
larger by a factor of 6.1'^\ ^^^ compared to the lower -density 



field at the same redshift. Digby-North et al. ( |2010 1 studied 
the z = 2.30 protocluster in the field of QSO HST700+643 
and found five protocluster AGN in their ^ 200ks Chandra 
observation. They identified AGN with Lx.h > 4.6 x 10"*^ 
erg/s in a sample of members selected via the BX/MD method 
and as LAEs. The X-ray AGN fractions in the BX/MD and 



7. SUMMARY 



LAE samples are 6.9'!^l and %, respectively, which 

are greater than the field fractions in similar galaxies at this 
redshift, particularly for the BX/MD sample. The protoclus- 



ter and field fractions mea sured by Lehmer et al. ( 2009) l and 
[Digby-No rth et al.| ( |2010 , for BX/MD galaxies) are shown 
on Figure [8[ While the field fractions for these high-redshift 
fields are measured in their immediate vicinity, this should 
be less important than for our sample because both less time 
has elapsed for environmental pre-processing and the environ- 
ments are not as evolved. 

These data on the evolution of the AGN fraction in the 
field, clusters, and protoclusters, while heterogeneous with 
respect to X-ray sensitivity thresholds, cluster selection, and 
host galaxy properties, broadly indicate that luminous AGN 
are anticorrelated with density in the local universe, found 
in cluster and field galaxies with approximately equal fre- 
quency at 1 < z < 1.5, and are correlated with local density 
by z > 2. This relative evolution of the fraction of lumi- 
nous AGN in field and cluster galaxies is consistent with the 
now-conventional picture that the most luminous AGN are fu- 
eled by gas-rich galaxy mergers ( Sanders et al.|1988 1 and the 
steady decline in the cold gas content of galaxies from high 
redshift to the present day. At the highest redshifts discussed 
here, which correspond to the redshifts of the protoclusters 
at2 < z < 3.1, a larger fraction of the baryonic mass fraction 
of massive g alaxies is in molecul ar gas compared to their lo- 
cal analogs ( Tacconi et al.[[20T0 i. The overdense protoclus- 
ter environment leads to greater likelihood of interactions and 
mergers of these gas-rich galaxies and thus the fueling of lu- 
minous AGN. By 1 < z < 1.5, the most overdense cluster en- 
vironments have grown substantially and are large enough to 
have substantial hot gas reservoirs, as illustrated by the detec- 
tion of extended X-ray emission from some clusters at z > 1 



( [Gobat et al.|20n]|Brodwin et al.|201 1| [Stanford et al.|2012| i. 
While there is still substantial star formation in many cluster 
members, including clusters in the present sample (Brodwin 
et al. 2013, in preparation), many other cluster galaxies have 
stopped substantial star formation, as indicated by the pres- 



ence of a color-magnitude relation (Mei et al. [2009 1. This 
transition is likely because they have largely exhausted their 
cold gas supply. The increase in the cross section for galaxy 
interactions and mergers in the cluster environment is conse- 
quently counterbalanced by the smaller fraction of gas-rich 
galaxies, e xcept perhaps on the outskirts of clusters (^Waggj 
[et al.|2012) l. Even at redshift 0.9 < z < 1.6 there is som e evr 



dence for an enhan cement of AGN in the infall region ( Pass 
bender et al.|2012 1. By the present day, cluster g alaxies have 



much less cold gas compared to field galaxies (Giovanardi' 
|et al.[[1983j [Oosterloo et al.,2010| ) and the relative velocities 
of their member galaxies are too great to produce bound pairs 
that eventually merge. Aside from the central cluster galaxy, 
which is likely fueled by gas cooling from the ICM, luminous 
AGN may only be found in pr esent-day clusters d ue to the in- 
fall of gas-rich field galaxies ( Haines et al.|2012| . This result 
is similar to the properties of star forming galaxies in clusters 



We have investigated the X-ray and MIR-selected AGN 
population in a sample of 13 clusters at 1 < z < 1.5 identi- 
fied by the ^p/fze r/IRAC Shallow Cluster Survey of [Eisen-] 
[hardt et al.| ( [2008[ ). We find a total of 1 1 X-ray counterparts to 
cluster members, eight of which have spectroscopic redshifts. 
There are also 27 MIR AGN associated with cluster members, 
12 of which have spectroscopic redshifts. All but three of the 
X-ray AGN are also MIR AGN. 

The X-ray AGN are quite luminous. All of the X-ray AGN 
have rest-frame, hard X-ray luminosities of Lx,h > 10*^ erg/s 
and four have Lx,h > lO^'* erg/s. AGN at these luminosities 
are extremely rare in low-redshift clusters, and in fact none 
have been found above Lx,h > 10'*'* erg/s in a combined sam- 
ple of44clustersatz< 0.5 that we constructed from studies 
by Martini et al. ( 2009| l and Haines et al. ( 2012| l. These new 
observations demonstrate that the order of magmtude increase 
in the cluster AGN fraction from the present to z ~ 1 continues 
toz-- 1.5. 

We have used photometric redshift estimates for galaxies 
out to five times the expected r20() radius of the clusters, or 
R = 2' - 10' from the cluster centers, to characterize the field 
population and likely field contamination due to the use of 
photometric redshifts to estimate cluster membership. These 
observations clearly indicate a substantial excess of galaxies. 
X-ray AGN, and MIR AGN associated with the clusters. We 
calculate the X-ray and MIR AGN fractions of all galaxies 
brighter than MJg+ 1, where MJf, corresponds to the knee of 
the luminosity function at the cluster redshift. We find that 
the field AGN fractions, defined from the sample of galaxies 
at the cluster redshift, but from R = 2'-6' (X-ray AGN) or 
R = 2' - 10' (IR AGN) from the cluster center, are consistent 
with the cluster AGN fractions within R < 2'. This stands in 
sharp contrast to estimates at low redshift, where the luminous 
X-ray AGN fraction is substantially lower in clusters relative 
to the field (e.g. Martini et al. 2009 ; Haggard et al. 201 0). 

The order of magnitude evolution of the cluster AGN frac- 
tion from z = to z ~ 1 .25 is greater than the rate of evolution 
of the field AGN fraction. While the luminous AGN fraction 
is approximately six times higher in the field than in clus- 
ters in the local universe, we find comparable fractions in the 
field and clusters at z ^ 1.25. Studies of two protoclusters and 
samples of field galaxies at even higher redshifts indicate that 
the luminous AGN fraction is higher in protoclusters than th e 
field at z > 2 ( [Lehmer et al.|2009[[Digby-North et al.|2010[ ). 
Taken together, these studies and our own demonstrate that 
there is a reversal of the local anticorrelation between lumi- 
nous AGN and local density at high redshift. The relative evo- 
lution of the AGN fraction in the field and clusters is strong 
evidence of environment-dependent AGN evolution. 
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